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, Using Monte Carlo simulations we study the dynamics of three-dimensional Ising models with 

fSJ . nearest-, next-nearest-, and four-spin (plaquette) interactions. During coarsening, such models de- 

» I ' velop growing energy barriers, which leads to very slow dynamics at low temperature. As already 

CIh, reported, the model with only the plaquette interaction exhibits some of the features characteristic 

of ordinary glasses: strong metastability of the supercooled liquid, a weak increase of the charac- 
teristic length under cooling, stretched-exponential relaxation and aging. The addition of two-spin 
interactions, in general, destroys such behaviour: the liquid phase loses metastability and the slow- 
dynamics regime terminates well below the melting transition, which is presumably related with a 
certain corner-rounding transition. However, for a particular choice of interaction constants, when 
. the ground state is strongly degenerate, our simulations suggest that the slow-dynamics regime ex- 

O ' tends up to the melting transition. The analysis of these models leads us to the conjecture that in 

the four-spin Ising model domain walls lose their tension at the glassy transition and that they are 
basically tensionless in the glassy phase. 
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I. INTRODUCTION 



A lot of efforts have been devoted in the last twenty years to understanding the behaviour of various glassy and 
disordered systems Such systems, which include conventional glasses, spin glasses, amorphous semiconductors, 
and many others are of great importance both experimental and theoretical. However, despite intensive research, our 
^ ' understanding of such systems is still limited. For example, even the very nature of the glassy phase in spin-glasses 
Q is still a very controversial issue Although they are much more abundant, conventional glasses seem to pose 

O , even a greater puzzle. Why do supercooled liquids fall out of equilibrium at a more or less well defined temperature? 
Why do they collapse into the glassy state when the cooling is fast enough and into the crystalline phase when the 
' cooling is slow? These fundamental questions still await definitive answers. One of the important problems in physics 
^ of conventional glasses is the continuing lack of a satisfactory microscopic model of such systems. In this respect the 
situation is much better for spin-glasses where it is commonly accepted that models containing quenched disorder 
, correctly describe physics of such systems. Lattice realizations of such models are a particularly valuable source of 
' information about spin-glasses The most realistic models of conventional glasses, so-called off-lattice models, still 
[ constitute an enormous computational challenge although progress in this field is also being made |^ . 

A model of conventional glasses should be capable of describing (at least) three phases: liquid, glass and crystal. 
I The actual state of the system should be determined by control parameter(s) (e.g., temperature) and possibly also 
"Y"; I its history. Since the glass is regarded as a liquid trapped during the falling out of equilibrium, the model should 

■ possess such a trapping mechanism. In spin glasses the trapping mechanism is related with energy barriers generated 
^ , by quenched disorder On the other hand we do not expect the quenched disorder to be a relevant factor in 

I ' conventional glasses because models with strong quenched disorder are unlikely to exhibit periodic solutions (which 
' O ] are needed for the model to be in the crystal phase). Recently, various lattice models, which do not contain quenched 
5 ' disorder, were studied which have some features of conventional glasses. Some of these models are infinite-dimensional 
and their thermodynamical properties can be found exactly |7|]. There are also finite dimensional models whose 
dynamics exhibit some glassy behaviour [^,^ 
. !^ Recently, it has been shown that the three dimensional Ising model with the four-spin (plaquette) interaction also 

^ . exhibits some glassy features |lo|-|l2|. This model undergoes a first-order phase transition between low-temperature 

■ (crystal) and high-temperature (liquid) phases. However, when conventional simulation techniques are used, the 
transition is screened by a very strong metastability during heating as well as cooling. For temperatures lower than 
the limit of metastability of the liquid phase, the model has a very slow coarsening dynamics. In addition, the 
zero-temperature characteristic length increases very slowly as a function of the inverse cooling rate, which is also an 
expected property of glasses. Further evidence of the glassy behaviour in this model has been recently reported by 
Swift et al. [n3| . They have shown that the glassy transition coincides with the divergence of a certain relaxation time 
and that aging properties of the model are also typical of glassy systems. They have also observed that some time 
correlation functions may decay as stretched exponentials. 

These results strongly suggest that the model with four-spin interactions might describe important aspects of the 
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glassy transition. It would be interesting to find which properties of this model are responsible for such a behaviour. 
It has been already suggested [p^|-Jl2| that the trapping mechanism might be related with diverging energy barriers. 
These barriers would arise in this model basically due to the same mechanism as in a model with competing nearest- 
neighbour and next-nearest-neighbour interactions examined by Shore et al. |^ (the SS model for short). However, 
the behaviour of the SS model is not fully consistent with our conception of glasses since it orders too quickly under 
cooling js). It was also suggested that the difference in the behaviour of the SS and four-spin models might be related 
with the degeneracy of the ground state in the four-spin model. This degeneracy might lead to some entropy barriers, 
which would be responsible for the strong metastability of the liquid phase. 

In the present paper, using Monte Carlo simulations, we examine a certain class of three-dimensional Ising models 
which generate energy barriers. These models are described by the following Hamiltonian 

H = — Ji ^ SiSj — J2 ^ SiSj — J4 ^ SiSjSkSi- (1) 

<ij> <<i,j» lid,k,l] 

In the above expression < .. > and << .. >> denote pairs of nearest and next-nearest neighbours, respectively, and 
[i,j,k,l] stands for summation over elementary plaquettes. In general, these models have double degenerate ground 
state and our simulations suggest that the dynamical properties in this case are similar to the SS model. However, 
when the interaction constants are such that the model has a strongly degenerate ground state (gonihedric case) , the 
dynamical properties change. Simulations suggest that two types of dynamical behaviour appear. In the first type the 
model behaves similarly to the already described four-spin model. In the second type, the glassy transition appears 
to coincide with the thermodynamic transition. 

Such behaviour gives rise to the following questions: why a glassy transition appears only in certain systems with 
slow dynamics and what is its nature. Analysis of the ground-state structure and thermodynamic properties of 
models studied here prompts the following conjecture, which, if confirmed, would constitute an important result of 
the present paper: at the glassy transition the domain walls lose their surface tension, and, as a result, the glassy 
phase is composed of tensionless domains. Although based on the analysis of Ising models, we hope that such an 
interpretation might shed some light on the nature of the glassy transition in more realistic systems also. Moreover, 
such an interpretation of glassy phase is in accord with some recent hypotheses concerning the nature of the glassy 
phase in spin glasses [Q. 

In section II we discuss briefly the properties of the four-spin model. In section HI we present the results of our 
simulations for the gonihedric case. The case of the doubly-degenerate ground state is discussed in section IV. Section 
V contains a summary of our results and some arguments on the nature of the glassy transition. 



II. ISING MODEL WITH PLAQUETTE INTERACTIONS 



This model corresponds to the case Ji = J2 = 0, J4 = 1, and has been already studied using Cluster Variational 



Method 1141 and Monte Carlo simulations 115 10 
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Clearly, the ferromagnetic configuration is a ground-state configuration of this model. It is also easy to realize 
that flipping coplanar spins does not change the energy. Thus any conflguration obtained from the ferromagnetic 
configuration by flipping coplanar spins is also a ground-state configuration. Moreover any combination of such 
coplanar flippings (even for crossing planes) does not increase the energy. Simple analysis along these lines shows that 
for the model on the lattice of the linear size L the degeneracy of the ground state is equal to 2^^^ . Although ground 
state of this model is strongly degenerate its ground-state entropy is zero. 

The model undergoes a first-order thermodynamic transition at T = Tc 3.6 which is, however, screened by very 
strong metastability As a result, when heated or cooled, the transition observed in simulations is shifted to 

T ^ 3.9 or T ~ 3.4, respectively. Transitions at these spinodals are accompanied by peaks in the specific heat. 

The low-temperature spinodal T 3.4 seems to coincide with the glassy transition. Below this temperature 
the model exhibits very slow coarsening dynamics p^ ] as well as aging properties which are characteristic of glassy 
systems [|l3| . A certain characteristic time, which governs the relaxation of energy-energy correlation functions, also 
seems to diverge at this temperature [ p3[ . In addition, the behaviour of the model under continuous cooling supports 
the glassy-transition interpretation of this temperature [C2I . 



III. GONIHEDRIC ISING MODEL 
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A. Ground state and thermodynamics 



It has ah'eady been suggested that the slow dynamics of the four-spin model might be related with energy barriers 
generated in that model These barriers arise due to the shape-dependence of the energy of excitations: it is 

not only the size of an excitation which determines its energy but also its shape. Such shape dependence appears also 
in the SS model. 

Are there any other models which could have a similar property? In our opinion, the shape dependence of energy of 
excitations should be rather a robust feature of Ising-type models. It is only in some specific cases, like the standard 
nearest-neighbour case, when this energy does not depend on the shape of an excitation. In particular, energy barriers 
appear in model (^) . The Hamiltonian of this model is quite general and it includes both the four-spin model and the 
Shore et al.'s model ( J4 = 0, Ji > 0, J2 <0). 

In the present section we examine a class of models described by this Hamiltonian, namely gonihedric models ||l6|,0 . 
These models correspond to the following choice of interaction constants: Ji — 2fc, J2 = and J4 = i(fc — 1). 
Gonihedric models have a strongly degenerate ground state. In addition to the ferromagnetic ground state any 
configuration obtained by flipping coplanar spins also minimizes the Hamiltonian. Any combination of such flips does 
not increases the energy, provided that flipping planes do not cross. As a particular example of such a ground-state we 
can mention lamellar configurations where e.g., every second plane of spins is flipped. Although lamellar structures 
constitute a legitimate ground-state, they do not survive at finite temperature as shown by Cirillo et al. using the 
Cluster Variational Method ||T^. We will return to this feature in the last section. 

For fc = the gonihedric model is equivalent to the four-spin Ising model. In this case the model has an additional 
symmetry which implies a larger degeneracy of the ground state since the flipping planes can now cross. As a result 
we obtain that antifcrromagnetic configurations belong to the ground state. Further analysis of differences between 
the fc = and the fc 7^ cases is postponed to the last section. 

Gonihedric models are expected to undergo a thermodynamic transition which for fc < fctr is of first order and for 
fc > fctr is of second order. Only very rough estimations of fctr(~ 0.5) are known p^ . In this section we analyze 
dynamical properties of the gonihedric model for k — 2, i.e., for a value with a continuous transition. Our results 
were obtained using a standard Monte Carlo method with random sequential update using Metropolis algorithm [ p^ . 
Some details can be found elsewhere ||ll|,|l^ . 

To find the thermodynamic transitions we measured the specific heat. Our simulations, which were made for various 
linear sizes L up to i = 40, locate the peak at Tc ~ 2.35, which is a good agreement with the Cluster Variational 
Method estimation jl^. The absence of hysteresis effects confirms that the transition at T = Tc is continuous. 
(Although the nature of the thermodynamic transition in gonihedric models is an interesting and still open problem, 
its further analysis is not an objective of the present paper). 



B. Dynamics 

An important indication of glassy dynamics is a slow evolution of a random quench. For usual models with 
nonconservative dynamics one expects pj| that the characteristic length / increases with time t as / ~ t^/^. However, 
in glassy systems I should increase much more slowly and presumably only logarithmically with time {I ^ \nt). Such 
a behaviour most likely appears in the SS model and in the four-spin model. In the following we present the results of 
our simulations of the evolution of quenches in the fc = 2 model. We measured the excess energy SE = E — E^o, where 
Eoo is the equilibrium energy. Our results for temperatures T < Tc are shown in Fig. ^ To relate the characteristic 
length with the energy excess we can employ the frequently used relation 

Ir^l/SE. (2) 

With this identification from Fig. ^ we infer that for all the examined temperatures the asymptotic increase of I is 
much slower than i^/^. At the end of this section we will argue that the relation (|^) most likely does not hold for this 
model and Fig. |^ actually suggest that the increase of / is slower than t^^^. Since there are no theoretical arguments 
for such a slow algebraic increase in our opinion it is quite plausible that asymptotically we have I ^ Int. Such a slow 
increase of / is most likely due to energy barriers. 

An alternative technique to examine the dynamics of the model is to measure characteristic times of certain 
processes. For example, one can measure the average time r needed for the inversion of a cubic like excitation. In 
this method, which was already applied to similar models PJll|, one prepares the system of the size L with fixed 
boundary conditions and interior spins which are opposite to the boundary spins. One expects that after some time, 
the system will invert the interior spins. For a two-spin Ising model or SS model above the corner rounding transition, 
T ~ L^, which indicates a relatively fast dynamics (naive inversion of this relation gives / ~ t^^^). On the other hand. 
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in the SS model below the corner rounding transition and in the four-spin model r increases much faster, presumably 
exponentially, with L. 

We measured the time needed for magnetization of the interior spins to reach the equilibrium value at a given 
temperature and the results are shown in Fig. |^. These results show that t even at the highest examined temperature 
increases faster than L^. This is a potential indication of an exponential increase t ~ a^{a > 1) in the entire 
low-tcmperature phase. 

Additional confirmation of such behaviour is obtained from simulations of this model under the continuous cooling. 
Similarly to simulations of the four-spin model |]l2| , we relax the random sample at a temperature Tq > and then 
continuously lower the temperature according to the formula T{t) = Tq — rt, where r is the cooling rate. When 
the temperature is reduced below the critical point Tc the growth of order begins. The slower the cooling the more 
ordered is the system at the end of the cooling , i.e., at T = (see Fig. |3|). To quantify the zero-temperature order 
we measure the excess energy SE at T = and the result is shown in Fig. ^ 

Using the relation (|^) this data suggests that asymptotically I ~ r~^^'^. Such a relation, which indicates that the 
growth of order is relatively fast, holds for the two-spin Ising model ||2^ and also for the SS model ||^. However, this 
conclusion is based on the validity of the relation and, similarly to the four-spin model we want to argue 
that this relation does not hold. Our argument refers to the following property of all gonihedric models: the energy 
of cubic-like excitations scales as their linear size L j2^. Let us recall that in two-spin Ising model, this energy scales 
as the area of the excitation (~ L^). Provided that the final configuration is composed of such domains of the size L 
and using a simple dimensional argument we obtain 

(3) 



A direct confirmation of the assumption about the structure of the configuration at the of the cooling process comes 
from visual inspection. In Fig. ^ we can see an example of single-layer configuration. One can clearly see cubic-like 
(i.e., non-rough) domains whose energy scales linearly with their size. Using the relation (||) Fig. ^ shows that the 
zero-tempcrature characteristic length scales as r^^/* which is much slower than in the two-spin Ising model but 
faster than in the four-spin model. 

In this section we used three independent techniques to probe the dynamics of the gonihedric model in the case 
of continuous thermodynamic transition k — 2. Domain coarsening suggest that the model has a slow dynamics up 
to, at least, the temperature T — 1.9 (Tc ^ 2.35). An analysis of the size dependence of the characteristic time r 
suggests that cubic-like domains remain non-rough at least up to the temperature T = 2.1 Thus, a slow-dynamics 
regime is most likely extended up this temperature. Since this is very close to the critical point is it not unlikely that 
a slow-dynamics regime actually covers the whole low-temperature phase. The behaviour of the model under cooling 
confirms such a scenario: if there would be a certain temperature Tq < T^ such that for Tq < T < T^ the dynamics 
would be fast then for the slow cooling the growth of order would be dominated by the time spent in this temperature 
interval and we would have I ~ r"^/^. Such a scenario takes place in the SS model The growth of order in our 
k — 2 model is much slower / ^ r^^/^ and excludes the existence of such a temperature Tq (unless it is very close 
to Tc and our simulations are not sufficient to detect the true asymptotic behaviour). Let us also notice that Shore 
et al. also analyzed certain SOS model for which (by necessity) Tq = Tc. Using some scaling arguments they have 
shown that for this model one should have I ^ r^^l^ ^ which is in agreement with our numerical result. 



IV. OFF-GONIHEDRIC ISING MODEL 

The gonihedric case corresponds to a certain choice of interaction constants in the Hamiltonian (Q). As we have 
noted this choice has important implications: the ground state is strongly degenerate and energy of excitations scale 
as their linear size and not as their area. In the present section we examine what is going on when the interaction 
constants of model (|l|) deviate from the gonihedric case p^ ]. 

As a particular example we choose: J\ =6, Ji — —1, J4 = 1/2, which differs from the gonihedric case fc = 2 by a 
modified nearest-neighbour coupling Ji . Such a model has a double-degenerate (ferromagnetic) ground state and our 
rough estimation of the critical temperature is Tc ~ 12.5. 

To examine the dynamics of this model we used the same techniques as described in the previous section. First, 
we examined the coarsening behaviour of this model. At low temperature (up to T ~ 3.0) we observed a very slow 
decrease of the energy toward the ground state value. Our data, which we do not present suggests that for such 
temperatures I most likely increases logarithmically with time. However, above this temperature dynamics becomes 
much faster and presumably the characteristic length increases as I ~ t^l"^ . 
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Such behaviour is confirmed by the measurements of the characteristic time r defined in the same way as in the 
previous section. Results of simulations are shown in Fig. ^. They indicate that for temperature T = 6 and 8 the 
characteristic time r increases as similarly to the SS model above the corner rounding transition. 

The above results indicate that the dynamical behaviour of the model is very similar to the SS model. Namely, in 
the low temperature regime the model has a slow dynamics and rapidly (faster than L^) increasing characteristic time 
T. However, within the ordered phase (i.e., for T < Tc ^ 12.5 there is also a high temperature regime where dynamics 
is much faster. Presumably, in this regime the dynamics is similar to other nonconservative systems with scalar order 
parameter . Similarity between the model examined in this section and the SS model is in our opinion related with 
the structure of the ground state: both models have only double degenerate (ferromagnetic) ground states. It is well 
known in the statistical mechanics that degeneracy of the ground state plays a very important role in determining 
the thermodynamic behaviour of models (e.g., critical behaviour). It is probable that this is also an important factor 
for determining the dynamical behaviour of models. 



V. NATURE OF THE GLASSY PHASE 

In the present paper we have examined the dynamics of models described by the Hamiltonian (^. Depending on 
the interaction constants we can distinguish three types of behaviour which are schematically shown in Fig. 0. 

(a) Double-degenerate ground state 

This, the most typical situation, appears in the off-gonihedric model studied in previous section and also in the SS 
model whose dynamics has been already examined in great details . The dynamics of the model a low temperatures 
(T < Tc) has two regimes separated by a certain temperature Tcr- For T < Tcr the model has slow dynamics 
with most likely logarithmically increasing characteristic length I. Such behaviour is related with the fact that 
at such temperature the model is below the corner-roughening transition p3|] . As a result, an evolving quench 
develops c omp licated structures of cubic-like excitations which are very stable and effectively block further coarsening 
dynamics |24| . At T = Tcr the model undergoes the corner rounding transition and the blocking mechanism is no 
longer effective. As a result the fast (standard) dynamics is restored. 

(b) Gonihedric case with continuous transition 

In this case the entire low-temperature phase has slow dynamics, whose origin is similar to the case (a). Namely, a 
quench develops cubic-like structures which block further coarsening dynamics. The degeneracy of the ground state 
seems to be the most important difference between this and the off-gonihedric case. Thus, we relate the disappearance 
of the corner-rounding transition (or maybe its overlap with Tc) with the infinite degeneracy of the ground state. 

(c) Four-spin model 

As we already mentioned, the four-spin model undergoes a dynamic transition which exhibits a lot of similarities with 
the glassy transition. It might be interesting to examine whether such a behaviour appears only at fc = or persists 
also for some other (small) values of k. 
Glassy transition: loss of surface tension 

Why does this transition exist in the four-spin case and not in the other cases? It was already suggested that the 
difference between the four-spin model and the SS model is related to the degeneracy of the ground state jl^. However, 
the above analysis of the gonihedric model with k — 2 shows that the infinite degeneracy of the ground state is not 
sufficient for the model to have a glassy transition (of course, we limit our analysis to models which can generate 
diverging energy barriers and thus have slow coarsening dynamics). Why does the gonihedric case k = 2 differ from 
the fc = case (i.e. the four-spin model)? Both models have strongly degenerate ground state. The degeneracy equals 
2'^^ for fc = and 2^ for fc = 2. Since in both cases degeneracy increase exponentially with the linear system size this 
difference does not seem important. In our opinion, however, the difference in the dynamical behaviour is related with 
the ground state structure of these models. As we already mentioned, for fc = the flipping planes, which generate 
various ground-state configurations might cross. As a result, in addition to ferromagnetic-like configurations we obtain 
antiferromagnetic-like ones. For the fc 7^ such crossings are not allowed and only ferromagnetic-like configurations 
are possible (we consider lamellar configurations also as ferromagnetic-like). 

This difference has important implications. Let us note that for the four-spin model in addition to tensionless 
domain walls, which appear for example when a cubic ferromagnetic 'up' domain is surrounded by 'down' one, there 
are tensionful ones too. As an example of such a domain wall we can consider an antifcrromagnctic domain surrounded 
by ferromagnetic one p^ . At first sight it does not seem to be much different from the k — 2 case. Indeed, when one 
considers a lamellar configuration where successive layers are of opposite sign (see Fig. ^) which is surrounded by a 
ferromagnetic domain than the excess energy of such a configuration scales as the area of the wall (i.e., the domain 
wall is tensionful). There is, however, an important question: do such configurations affect coarsening dynamics or, 
in other words, are they spontaneously generated in sufficient amounts? In our opinion the answer to this question is 
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negative. Our first argument comes simply from tfie visual inspection of the snapshot configuration. In Fig. ^one can 
see relatively large ferromagnetic-like domains but there is no indication of lamellar ones. The second argument comes 
from the Cluster Variational Method calculations by Cirillo et al. who have shown that the lamellar structures are 
equivalent (i.e., of the same energy) to the ferromagnetic one but only for the ground state. At non-zero temperatures 
they are always metastable. These arguments show why such configurations are not spontaneously generated during 
the evolution of the quench. They imply that for k = 2 the dominant domain walls which exist at the late-time 
evolution of the quench arc tensionless. 

On the other hand, for the four-spin model, antiferromagnetic structures are fully equivalent to the ferromagnetic 
ones. In the liquid phase both ferromagnetic and antiferromagnetic domains are intertwined and form very complicated 
structures. As already noticed fl^ , tensionful domains usually have lowest energy barriers and the system can 
relatively easily remove the interior domains. On the contrary, tensionless domain walls have large energy barriers 
and their dynamics is much slower. It means that in the liquid phase dynamics is dominated by dynamics of tensionful 
domains and thus resembles the dynamics of two-spin Ising models. 

It is in our opinion very likely that upon lowering the temperature the system will undergo a phase transition which 
will eliminate tensionful domain walls. Below that transition the energy of the system would be located mainly in 
tensionless domain walls. It means that at this transition the antiferromagnetic-ferromagnetic symmetry of the model 
would be spontaneously broken. In other words, at this transition the system selects a dominant type of domains, 
whether ferromagnetic or antiferromagnetic. This, rather novel, type of symmetry breaking is an essential ingredient 
of the transition which we tentatively identify as the glassy transition. 

Finally, let us note that the idea that the glassy phase consists of a complicated mixture of tensionless domain walls 
appeared recently in the context of spin glasses |Q. It suggests that, at least at the geometrical level, spin glasses and 
ordinary glasses might exhibit a lot of similarities. Their further explorations is, however, left as a future problem. 
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FIG. 1. The excess energy 5E as a function of t in the log-log scale. Simulations were made for the system size L = 100 
and T = 1.0, 1.3, 1.6 and 1.9 (from bottom to the top). The dotted line has a slope 0.5. 

FIG. 2. The characterisitc time r needed for the inversion of the cubic excitation of the size L as a function of L for the 
gonihedric case (A; = 2). Simulations were made for T = 1.9 (o), 2.1 (□) and 2.3 (V). The dotted line has a slope corresponding 
to r ~ L^. 



FIG. 3. The energy _B as a function of temperature for (from the top) r = 0.02, 0.002, 00002 and 0.00002. The ground-state 
energy equals (for k = 2) -4.5. 



FIG. 4. The excess energy SE as a function of r in the log-log scale. The dotted line has slope 0.45. 

FIG. 5. A single-layer snapshot of a Monte Carlo configuration at the end of the cooling process. Simulations were made 
for L = 50 and r = 0.001. 

FIG. 6. The characteristic time r needed for the inversion of the cubic excitation of the size L as a function of L for the 
off-gonihedric case. Simulations were made for T = 2 (0), 3 (A), 6 (□) and 8 (o). The dotted line has a slope corresponding 
to r ~ L^. 

FIG. 7. Three types of dynamical behaviour found for model (0). (a) Double-degenerate case (SS model); (b) the gonihedric 
model with continuous thermodynamic phase transition; (c) the four-spin model (and possibly the gonihedric model with 
discontinuous thermodynamic phase transition). 

FIG. 8. A two dimensinal (4x4) lamellar domain surrounded by the ferromagnetic domain. The excess energy comes only 
from horizontal boundaries of domain wall and thus scales as the size of the interior domain. In the three-dimensional case the 
excess energy would scale as an area of interior domain. 



7 








0.0 



0.5 1.0 

iogio(L) 



1.5 



Fig.2 




-5 -4 -3 -2 -1 



login(r) 



Fig.4 



!»>:»::> 
:::::::::: 



88888888888:::::::::::::::::::::::! 



!:::::::::^ 



;:::::::::::::::! 



!:::::::88^ 



::::::::8888:::::::::::::::i 



:::::::::::::::::::::^ 



!:::::::s 



until 



::::::s 



!:::::::s 



;88888:::ssssss::::::::::::::i 



!:::::::::^ 



10 



!8:::::::::88:::::::::::::::::88888:i 



20 30 

lattice sites 



40 



50 



Fig.5 













O 


v 














o 


o 




A 


A 






o 


o 

A 

' 8 


o 

A 
O 


A 

o 


A 

□ 

o 


□ 
O 


□ 
O 




o. ■ ■ ' 
,.■■6 


















1 1 1 



0.0 0.5 1.0 1.5 2.0 



iogio(L) 



Fig.6 



Tcr Tc J 

n 1 \ ^ 

(a) 

Tc(=Tcr) T 

-n 1 ^ 

(b) 

Tg Tc T 

1 ^ \ ^ 

(c) 



Fig.7 




Fig.8 



